Self-organized periodicity of protein clusters in growing bacteria 
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Chemotaxis receptors in E. coli form clusters at the cell poles and also laterally along the cell 
body, and this clustering plays an important role in signal transduction. Recently, experiments 
using fluorescence imaging have shown that, during cell growth, lateral clusters form at positions 
approximately periodically spaced along the cell body. In this paper, we demonstrate within a 
lattice model that such spatial organization could arise spontaneously from a stochastic nucleation 
mechanism. The same mechanism may explain the recent observation of periodic aggregates of 
misfolded proteins in E. coli. 
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Spatial organization of proteins is important in many 
cellular processes including growth, division, movement, 
and establishment of polarity [1]. Over the past few 
years, advances in imaging techniques such as fluores- 
cence microscopy have led to an increased appreciation 
of the scope and character of protein organization in 
cells. For example, in Escherichia coli and other bac- 
teria, chemosensory complexes form large clusters con- 
taining thousands of receptors [2]. Clustering of these 
receptors plays a crucial role in the signal integration 
and receptor cooperativity required for chemotaxis [3], 
i.e. directed movement in chemical gradients. Recent 
work by Thiem et al. |H [5] demonstrated that clusters 
of chemotaxis receptors are approximately periodically 
positioned along the cell wall, independent of any known 
positioning mechanism such as the Min system jB] . Other 
examples of periodically positioned protein clusters have 
emerged as well [311]; for example, protein aggregates as- 
sociated with cellular aging in bacteria exhibit a regular 
distribution along the the cell's long axis in filamentous 
E. coli 8 . The question arises - could such periodic 
positioning arise spontaneously or does it require the ex- 
istence of an unknown positioning system? 

Here we demonstrate, within the context of a minimal 
lattice model, that protein clustering and periodic po- 
sitioning of clusters can emerge spontaneously in grow- 
ing cells. Lattice models have been used before to study 
clustering of membrane proteins with short-range inter- 
actions [9]. In our model, existing clusters act as sinks 
for proteins newly inserted in the membrane, locally re- 
ducing the density of protomers and thus preventing nu- 
cleation of new clusters. As cells grow, existing clus- 
ters separate, ultimately allowing new clusters to nucle- 
ate at a characteristic spatial separation set by insertion, 
diffusion, interaction strength, and growth rates. The 
proposed mechanism is quite general; while we focus on 
membrane proteins, the mechanism also applies to ag- 
gregation of cytoplasmic proteins in the body of the cell 
(e.g. misfolded protein aggregates) [3]. 

In our model, the cell membrane is represented by a 



square lattice, whose a;-axis coincides with the long axis 
of the cell (see Fig. [TJ. We employ periodic boundary 
conditions in the y direction to account for the cylindri- 
cal shape of bacteria like E. coli. The protomers (inde- 
pendently diffusing protein units) associated with the cell 
membrane |10j are treated as particles which can perform 
random walks on the lattice. Each lattice site is therefore 
associated with a variable cr^, either occupied, cr.j = 1, or 
empty, o~i — 0. We assume a nearest-neighbor attractive 
interaction between particles with interaction energy J, 
measured in units of the thermal energy k^T. To control 
the nucleation barrier, we also include a conformational 
energy cost given by a J for each particle with any neigh- 
bors, which accounts for the loss of internal entropy when 
a particle associates with a cluster or a second protomer. 
The total energy of the system (in units of k^T) is 

E = — J o-iO-j + aJn c , (1) 

<i,3> 

where n c is the total number of particles with one or 
more neighbors (i.e., the number of particles that are 
in clusters of size two or greater) . Experiments indicate 
that the lateral receptor clusters are relatively immobile 
while individual membrane proteins are typically free to 
diffuse |4J. In our lattice model, we therefore consider 
only movement of individual particles. 

We use a Metropolis Monte Carlo algorithm to simu- 
late the system. A randomly selected particle is moved 
to one of its unoccupied neighboring sites with an ac- 
ceptance probability p = min(l, e~ AE ) where AE is the 
energy change due to the proposed displacement of the 
particle. One Monte Carlo time step corresponds to one 
attempted move for each particle present. 

In the absence of the conformational energy cost 
(a = 0), the thermodynamic system described by our 
energy function can be mapped to a two-dimensional 
Ising model, for which the critical interaction strength 
is known to be J c « 1.763 [11 . When the interaction 
strength is low, J < J c , the system has one stable homo- 
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geneous phase, while for J > J c , the system can phase 
separate into regions of high and low density. The con- 
formational energy cost increases J c . 

To account for growth of the bacterial cell, we allow 
the lattice to expand in the x direction according to 
L x (t) w Lj;(0)e 7t , where L x (0) is the initial length of 
the bacterium and 7 is the growth rate. The expan- 
sion of the lattice is implemented by random insertion of 
empty columns at a rate jL x (t) with equal probability 
anywhere in the lattice. Based on the observation that 
newly synthesized chemotaxis receptors are inserted into 
the cell membrane uniformly over the entire length of the 
cell JT2], particles are randomly deposited onto the lattice 
at a rate k on per available (i.e. unoccupied) site. This 
ultimately leads to an average density of occupied sites 
Po = kon/(l + k on ). To see this, let N(t) and n(t) be 
the total number of lattice sites and the number of occu- 
pied sites, respectively, at time t, with N(t) — iV(0)e 7t . 
On average, the total rate of particle deposition is given 
by dn/dt = k on [N(t) — n(t)]; the general solution to this 
equation is given by n(t) = p a N(t) + Ce~ kollt , where 
Po = fcon/(7+ fcon) is the asymptotic value of the particle 
number density, defined by p(t) = n(t)/N(t), and C is a 
constant. We start our simulations with p(t — 0) = p Q so 
that, on average, p(t) remains fixed at p a . 

For our simulations, the cell circumference was fixed 
at L y = 50, with a = 0.5 and interaction strength 
J = 4, considerably greater than the critical strength 
J c . The system was initialized with a cluster at each end 
of the cell to mimic the existing clusters at the poles. 
At the start of each simulation, the length of the cell is 
L x (0) = 20 and, as the cell grows, newly inserted par- 
ticles aggregate to form clusters that grow with time. 
Fig. 2 illustrates a series of snapshots from a representa- 
tive run with the growth rate chosen to be 7 = 0.8 x 10~ D , 
and a deposition rate k on = 2x 10~ 6 , yielding p a = 0.2. 
Notice that clusters of particles spontaneously appear 
at positions approximately periodically spaced along the 
cell. 

The emergence of a self-organized periodicity of clus- 
ters can be understood by noting the positions of new 
clusters when they first appear. At the start of the sim- 
ulation the clusters at each end of the cell act as sinks 
for newly inserted particles. As the cell grows and these 
two clusters move apart, a new cluster forms roughly at 
the midpoint of the cell. As cell growth continues, newly 
inserted particles spontaneously aggregate to form new 
clusters in between existing clusters; the location of a 
newly formed cluster is preferentially at the middle of 
two existing clusters, resulting in periodically positioned 
clusters (Fig. [2]). If the separation between two existing 
clusters is below a characteristic length £ c , diffusion dom- 
inates and particles are absorbed by old clusters; if the 
separation is larger than £ c , particles nucleate to form a 
new cluster, leading to periodic spacing ~ £ c . 

To quantify the separation between clusters, we ob- 
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FIG. 1: Schematic of the lattice model. Particles hop at 
random between neighboring lattice points and can join or 
leave an existing cluster from the boundary of the cluster. 
Columns of lattice points are inserted at random to mimic 
cell growth, and particles are inserted at random to mimic 
protein insertion in the membrane. 
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FIG. 2: Snapshots of the model cell membrane at times 70, 
75, 80, 85, and 90 x 10 4 Monte Carlo time-steps as defined 
in the text. The total size of the system in the final snap- 
shot is 1350 columns by 50 rows, with w 20% of lattice sites 
occupied by protein particles (green). In the inset we show 
a representative image of protein receptor clusters in E. coll 
from Victor Sourjik's lab (courtesy: Victor Sourjik) 



taincd the distribution of the separations between neigh- 
boring clusters for systems grown to L x — 1900 (see 
Fig. [3]). Due to stochastic fluctuations, a cluster is not 
thermodynamically stable until it reaches a critical size. 
We found that clusters of size > 50 were stable and did 
not disappear; we thus used size 50 as a criterion to iden- 
tify a cluster. The separation between neighboring clus- 
ters is defined to be the distance between the centers of 
mass of the two clusters. For L x = 1900, there are on 
average eight clusters present in the system. The dis- 
tribution of separations exhibits a single maximum at 
Ax ~ 230, indicating a preferred separation between 
neighboring clusters. Moreover, the fraction of cluster 
separations less than half the peak value is only 7.6%, 
indicating a strong suppression of close (i.e. Ax < 115) 
clusters. For comparison, the distribution of inter-cluster 
separation would be an exponential if the cluster centers 
were positioned randomly (dotted line in Fig. [3]). 

To investigate the mechanism responsible for cluster 
positioning, we studied how the position of a newly 
formed cluster depends on the positions of existing clus- 
ters. For simplicity, we used periodic boundary condi- 
tions in the x direction and initialized simulations with 
L x (0) = 1, that is, with only one column of sites in the 
system. We used a = 0.5, J = 4.0, k on — 1.1 x 10~ 6 , and 
7 = 10~ 5 , yielding po w 0.1. As the system grows, the de- 
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FIG. 3: The distribution of inter-cluster separations. Simu- 
lations were stopped when the system reached L x — 1900, 
which corresponds to approximately eight clusters in the 
system. The density is p — 0.1 and the growth rate is 
7 = lx 10~ 5 . The data was averaged over 70,000 simulations. 
The bin size is 20. If the centers of clusters were randomly 
distributed, the distribution of inter-cluster separations would 
be exponential as shown by the dotted curve. 
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posited particles aggregate to form first one cluster and, 
later, a second stable cluster. We recorded the position 
of the second cluster (once it had reached a size > 50) 
with respect to the first. The distribution of the separa- 
tions between the new cluster and the existing cluster as 
a fraction of the total cell length is plotted in Fig. 4(a). 
The distribution peaks at a reduced distance of 0.5, i.e., 
the second cluster is most likely to form at the midpoint, 
equidistant from the edges of the existing cluster. 

To understand why the second cluster formed near 
midcell, we investigated the density profile of particles in 
the dilute region between two stable clusters during cell 
growth. Simulations were performed, as for Fig. 4(a), us- 
ing periodic boundary conditions in x and starting from 
a single column. In Fig. 4(b), we plotted the particle 
density profile along the x axis as measured when the 
system reaches L x — 300 at an average global particle 
density p a = 0.1. In these simulations, at L x = 300, typ- 
ically there was one stable cluster in the system. Position 
was measured from one edge of the cluster and was nor- 
malized by the distance between the two flanking edges 
of the cluster. The average particle density fits very well 
to a quadratic function with the maximum at midcell. 

The observed quadratic particle density profile can be 
understood as follows. The average local particle density 
p(r,t) in regions that do not contain clusters satisfies 
the diffusion equation dp/dt — DV 2 p + k on , where D is 
the particle diffusion coefficient and k on is the particle 
insertion rate. Consider a region flanked by two stable 
neighboring clusters with flanking edges at x = and 



FIG. 4: (a) Distribution of the separation along the x di- 
rection between the center of the new (second) cluster and 
the center of the old (first) cluster, using periodic bound- 
ary conditions in x (see inset). The dotted line is a guide 
to the eye. (b) The particle density profile in the x direc- 
tion (periodic boundary conds.) with one cluster in the sys- 
tem. The particle density was measured when the system 
reached L x = 300 and was averaged over 10,000 simulations. 
The distance was normalized by the separation between the 
two flanking edges of the cluster as measured in the x direc- 
tion. The smooth line is a fit to a parabola, (c) The average 
inter-cluster separation versus the particle density p at fixed 
growth rate 7 = 10 -5 . (d) The average inter-cluster separa- 
tion versus growth rate 7 at fixed particle density p a = 0.1. 
The smooth curves in (c) and (d) are fits to power laws Po 47 
and 7 _0 - 53 respectively, (e) The cluster size distribution av- 
eraged over 10,000 simulations. Cluster sizes were measured 
when the system reached L x — 1900; the red line is a fit to 
the curve y = exp[aia; + a2a;lnx-|-a3 In x + 04], with ai « 0.8, 
a,2 ~ —0.37, as ~ —7.6, and 04 w —0.9, for cluster size, 
(f) Dependence of the standard deviation of cluster separa- 
tion (normalized by the average separation) on interaction 
strength J. The parameters and the method for collecting 
data were the same as in (a). In the inset we show the depen- 
dence of inter-cluster separation on J. 



x = I. Approximating stable clusters as perfect sinks for 
particles gives the boundary conditions p(0) = p{t) = 0. 
Hence, in the membrane strip between two clusters, the 
steady-state solution 13J to the diffusion equation is 

"« = -5iH) 2 + ¥- < 2 > 



4 



The peak of p(x) is located at x = 1/2, precisely at the 
midpoint between two cluster edges. The maximum par- 
ticle density is p max = k on £ 2 /8D. Notice that it scales 
quadratically with cluster spacing. 

Since the particle density peaks at the midpoint be- 
tween two neighboring clusters, this is the most likely 
place for a new cluster to nucleate. Of course, nucleation 
of a new cluster is a stochastic event; nevertheless, the 
probability of nucleating a new cluster is a highly nonlin- 
ear function of the local density, so the quadratic density 
profile gives rise to the sharply peaked distribution for 
the position of the new stable cluster (Fig. 4(a)). 

The existence of a relatively sharp density threshold 
for cluster nucleation predicts scaling relations for the 
mean cluster separation. For steady-state growth of a 
long cell, on average a new cluster must appear be- 
tween neighboring old clusters every time the cell doubles 
(e.g., see Fig. |2j. If there is a sharp density thresh- 
old pthrcsh for cluster nucleation, then a new cluster 
will nucleate when the peak density between old clus- 
ters reaches approximately this value. We can therefore 
estimate the upper limit of cluster separation in terms 
of Pthrcsh (with the lower limit being a factor of two 
smaller, and the mean lying between the two.) According 
to Eq. [2j the peak density between clusters depends on 
their separation £ according to p max = k on £ 2 /8D, where 
k on = p j/(l — Po) ~ p j. Nucleation of a new cluster 
should therefore occur when p max w Pthresh j that is, for 
a cluster separation £ c w 2[2Z?p t hrcsh/(Po7)] 1 ^ 2 - This im- 
plies that the mean cluster separation in a growing cell 
will obey the scaling relations £ c ~ p l ' 2 and £ c ~ 
[13] , which are verified in Figs. 4(c) and 4(d). 

We have demonstrated that an interplay of protein 
creation or insertion in the membrane, protomer diffu- 
sion, aggregation, and cell growth can lead to periodically 
spaced protein clusters. In addition to the larger stable 
clusters, our stochastic nucleation mechanism also gen- 
erates a quasi-steady state distribution of smaller clus- 
ters as shown in Fig. 4(e). Finally, in Fig. 4(f), we plot 
the dependence of the standard deviation of inter-cluster 
separation on the interaction strength J; the data implies 
that the periodic placement of protein clusters is robust 
to variations in J provided J < 5. 

For chemoreceptors in E. coli, the observed inter- 
cluster spacing is of the order of 1/xm. Assuming a 
membrane diffusion constant of D = 0.018pm 2 /sec [TS] 
and doubling time of around 60 minutes, we estimate 
Pthresh ~ 2.4 x lCP 3 p , implying that the vast majority 
of the chemoreceptors are bound to clusters rather than 
existing as free protomers. This prediction is testable us- 
ing modern fluorescent imaging techniques [16) . In prin- 
ciple, the scaling relations for £ c as a function of recep- 
tor density and growth rate can also be tested experi- 
mentally by measuring cluster spacings in cells over- or 
under-producing chemoreceptors and growing at differ- 
ent rates. Moreover, our model predicts that if proteins 



are produced (and hence inserted into the membrane) in 
bursts rather than continuously, it would adversely affect 
the periodic placement of clusters. 

Finally, we consider the biological significance of the 
periodicity of protein clusters. In the case of chemorecep- 
tors, periodic cluster spacing ensures that each daughter 
cell receives at least one cluster following cell division or 
fragmentation of a filamentous cell [I]. This "equiparti- 
tion" of protein clusters among daughter cells might be 
advantageous as well for other protein complexes. In con- 
trast, the larger spacing between misfolded protein ag- 
gregates typically means that only one daughter receives 
an aggregate; this unequal partition results in different 
fitnesses of the two daughters, with a potential overall 
reproductive advantage for the population [8]. 
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